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Abstract 

We study the nature of the synchronization transition in spatially extended systems by discussing 
a simple stochastic model. An analytic argument is put forward showing that, in the limit of 
discontinuous processes, the transition belongs to the directed percolation (DP) universality class. 
The analysis is complemented by a detailed investigation of the dependence of the first passage 
time for the amplitude of the difference field on the adopted threshold. We find the existence of 
a critical threshold separating the regime controlled by linear mechanisms from that controlled by 
collective phenomena. As a result of this analysis, we conclude that the synchronization transition 
belongs to the DP class also in continuous models. The conclusions are supported by numerical 
checks on coupled map lattices too. 
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I. INTRODUCTION 



Synchronization in dynamical systems has recently become the subject of an intensive 
research activity for various reasons that range from the application to transmission of 
information, to the spontaenous onset of coherent behaviour and also because because it is 
one of the mechanisms controlling the degree of order present in a chaotic evolution. Most 
of the attention has been, so far, focused on the behaviour of low-dimensional systems. As 
a result of these investigations, several kinds of synchronizations have been identified (the 
most important being phase and complete synchronization) and the corresponding transition 
scenarios characterized 

More recently, the interest has shifted towards high-dimensional chaos and, specifically, 
on the behaviour of extended systems, a context in which an overall picture is still lacking. 
In this paper, we devote our interest to complete synchronization in lattice systems. This 
kind of synchronization has been introduced and studied into two different setups. In the 
former one, identical copies of a given system (with different initial internal states) converge 
to the same trajectory, when forced with the same random signal. This, so-called stochastic, 
synchronization can occur only if the dynamics resulting from the stochastic forcing becomes 
linearly stable, i.e. the maximum Lyapunov exponent is negative 0, ^ ^, ^, 0. In the 
latter setup, two identical systems are coupled together: if the coupling strength is strong 
enough, both eventually follow the same, chaotic, trajectory. This is the so called chaotic 
synchronization. For it to be observed, it is sufficient that the transverse Lyapunov exponent 
is negative Therefore, in low-dimensional systems, synchronization transition can always 
be reduced to a linear stability problem. 

On the other hand, recent numerical investigations [|, ^ indicate that the synchronization 
scenario in spatially extended dynamical systems exhibits more complex and interesting 
features. In fact, the addition of the spatial structure may turn the linear stability problem 
of low-dimensional systems into a nonequilibrium phase-transition problem. 

In analogy to low-dimensional systems, various coupling schemes have been already con- 
sidered. For instance, "stochastic synchronization" has been studied in Coupled Map Lat- 
tices (CML's), by adding the same spatio-temporal noise ^{x,t) to different trajectories, 
Ui{x,t) and U2{x,t), of the same system p[, according to the following scheme 

Ui{x,t + 1)= f[ui{x,t) + VSi^,t)]+a^{x,t) , z = l,2 (1) 
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where 

V'lu{x,t) = -u{x + l,t) + -u{x — l,t) — eu{x,t) , (2) 

is the short-hand notation for the discretized Laplacian operator {e plays the role of a 
diffusion constant) and f[x] is a map of the unit interval able to generate chaotic behaviour. 
Moreover, a is the amplitude of the forcing term, x is an integer index labeling the lattice 
sites, t is a discrete time variable and the noise term is assumed to be bounded and 6- 
correlated in space and time, i.e. {^{x,t)C,{y, s)) oc S^^ySt^s- Synchronization is possible 
when the difference w{x,t) = \ui{x,t) — U2{x,t)\ between simultaneous configurations of 
the two systems converges everywhere to zero. The stability coefficient of the solution 
w{x,t) = is usually called the Transverse Lyapunov Exponent (TLE). In the context of 
stochastic synchronization, the evolution of a small w{x,t) reduces to the tangent dynamics 
of the single CML, so that the TLE coincides with the maximum Lyapunov exponent of the 
noise-affected dynamics. Accordingly, synchronization can arise only when the stochastic 
forcing induces a negative maximum Lyapunov exponent. This is possible if the probability 
distribution of the state variable mostly concentrates in the region of the interval where the 
map acts as a contraction. 

Alternatively, one can study the behavior of two directly coupled systems 0| 

ui{x,t + l) = {1 - a)f[ui{x,t) + Vlui{x,t)] + af[u2{x,t) + V%{x,t)], 

U2{x,t + 1) = {1 - a)f[u2{x,t) +V^^U2{x,t)] + (rf[ui{x,t) + Vlui{x,t)] . (3) 

At variance with the previous case, the coupling strength a modifies the evolution law of 
w{x,t), by adding a stabilizing term, while it leaves unaffected the dynamics of the fully 
synchronized regime. Accordingly, the TLE may become negative, while the maximum 
Lyapunov exponent, unchanged, remains positive. 

While the negativity of the TLE is always a necessary condition to observe synchroniza- 
tion in spatially extended systems, for smooth enough dynamical systems, it proves to be 
sufficient too. In fact, the study of stochastic and chaotic synchronization, carried on in 
Refs. Pi and |^, respectively, have shown that synchronization occurs as soon as the TLE 
becomes negative and, correspondingly, the propagation velocity of finite-amphtude per- 
turbation vanishes. In particular, Ahlers and Pikovsky argue that the dynamics of the 
coarse-grained absolute value w of the difference field is described by the following stochastic 
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partial differential equation, 

dtw{x,t) = DV'^w{x,t) + ciw{x,t) — C3W^{x,t) + w{x,t)r]{x,t) , (4) 

where D > 0, C3 > and the Gaussian noise term rj is (5-correlated in space and time, 
i.e. {ri{x,t)ri{y, s)) oc 6{x — y)6(t — s) . This equation is formally equivalent to the mean- 
field equation of the class of Multiplicative Noise (MN) nonequilibrium critical phenomena 
||T0| . By a Hopf-Cole transformation, h{x,t) = —lnw{x,t), the above equation can be 
transformed into |Tl[ 



dth{x, t) = DV^h{x, t) - D{Vh{x, t)f - (ci - ^) + cse-^'^^"'*) + 7]{x, t) (5) 

describing the critical behaviour associated with the depinning transition of a Kardar-Parisi- 
Zhang (KPZ) interface from a hard substrate. Numerical analysis confirms that the critical 
exponents evaluated for the two different coupling schemes are both compatible with those 
predicted for the MN model. 

On the other hand, it has been observed that in the presence of strong and localized non- 
linearities, the non- synchronized regime may coexist with a negative TLE [|], . In this case, 
the transition does occur when the propagation velocity of finite-amplitude perturbations 
vanishes, while its critical properties turn out to belong to the class of Directed Percola- 
tion. Such an equivalence has been found by noticing that the fraction of non- synchronized 
sites (defined as those points where is larger than some small fixed threshold) is the 

appropriate order parameter corresponding to the fraction of active sites in DP. 

In this case, one cannot follow the same derivation as above, because even close to the 
critical point, the evolution equation for w{x,t) cannot be linearized, since it is precisely 
the nonlinear effects which guarantee a propagation of finite-amplitude perturbations in the 
presence of a negative TLE. It is worth recalling that in the formulation of Reggeon Field 



Theory, the DP transition is described by the effective equation [T^, |T4 



dtp{x, t) = DV'^p{x, t) + cip{x, t) - C2p'^{x, t) + a/p(x, t)?7(x, t) , (6) 

where p{x,t) is the density of active sites and C2 > 0. Behind the similarity between this 
and Eq. (^), one should notice the crucial difference in the noise amplitude: the square-root 
versus linear dependence on p is indeed responsible for turning the MN critical behaviour 
into a DP-like one. In this paper, we plan to explain why the presence of a discontinuity 



(or a strong nonlinearity) may lead to the effective equation (^. To this aim, in Section 2 
we introduce a simple Random Multiplier (RM) model as an effective equation for the time 
evolution of the difference variable w{x,t) for discontinuous and strongly nonlinear CML's. 
This model was originally introduced in [|T3| to account for the mechanism of propagation of 
information in stable chaotic systems. We analyze its phase diagram, and we also discuss how 
the synchronization transition may be modified when a true discontinuity in the dynamics is 
changed into a strongly nonlinear continuous mapping. The relation between the RM model 
and the DP mean-field equation (^) is analyzed in Section 3. 

There is a further basic question that will be addressed here. All microscopic models that 
are known to exhibit a DP critical behaviour are defined by referring to discrete and finite 
state variables, such as the probabilistic cellular automaton model proposed by Domany and 
Kinzel In such cases, the so-called "absorbing state" can be unambiguously identified. 
For instance, in the cellular automaton of Ref. |T^, a sequence of "0"s can only change 
from its boundaries (this is the reason they are defined as contact processes). In the context 
of synchronization, the dynamical variable is continuous and the condition w{x,t) = is 
never exactly fulfilled at any finite time, even in a system of finite size. As a consequence, 
in numerical experiments 0, ^ one has to fix a small, but somehow arbitrary, threshold 
value, below which the trajectories are assumed to be synchronized. The same numerical 
simulations show that, independently of the dynamical rule, when the space average of 
w{x,t) decreases below a threshold value O{10^^) it does not grow again. However, one 
cannot a priori exclude that a large fluctuation of some local multiplier drives the system 
out of this weakly absorbing state. On the contrary, it looks plausible to assume that in 
an infinite system such a large fluctuation occurs with probability one. In Section 4 we 
tackle the problem of the existence of an effective absorbing state even in the presence of a 
continuous state variable. The study of the first passage time required for the space average 
of the difference variable w{x, t) to go through a series of decreasing thresholds clarifies that, 
contrary to intuition, it is possible to assign an effective finite "measure" to the synchronized, 
i.e. absorbing, state. Finally, conclusions are drawn in Section 5. 



5 



II. THE RANDOM MULTIPLIER MODEL: DEFINITION AND PHASE DIA- 
GRAM 



In this section, we introduce the RM model with the aim of closely reproducing the 
synchronization transition occurring in coupled piecewise linear maps of the type, 



xjai < X < «! 

1 — (a; — ai)/a2 ai < x < ai + a2 ■, (7) 
{x — ai — a2)l(X\ «! + a2 < x < 1 



where < ai < 1 and < 0:2 < 1 — «i. For any 0:2 > 0, the map is continuous with a 
highly expanding middle branch (when 02 ^ !)• In the limit 0:2 = 0, j[x) reduces to the 
discontinuous Bernoulli map with expansion factor \ja\. 

In the bidirectional synchronization setup (|^), the corresponding TLE is 

A± = AM + ln(l-2(T). (8) 

where Am is the maximum Lyapunov exponent of the single, uncoupled, chain. Therefore, 
linear stability analysis indicates that a small deviation w(a;, t) = \u\(x^t^ — ui[x^t^\ is 
contracted when a > (1 — 1/Aa/)/2. However, this is not the whole story even in the 
absence of multiplier fluctuations, because whenever t) and U2{x,t) fall on different 

sides of the map-discontinuity, w{x,t) becomes at once of order 1, being amplified by a 
factor close to l/w{x,t). The probability of such events depends on the probability density 
of the variables in the case of a sufficiently smooth distribution across the discontinuity. 



the probability is, to a leading order, proportional to w itself |T^. The same qualitative 
behaviour does occur also for «2 > 0, except that now, when w < a2, the amplification 
factor cannot be larger than (1 — 2a)a2- Moreover, the probability of such amplifications 
does no longer depend on w. 

In the following, instead of determining the local dynamics of w from the actual evolution 
of Ui{x,t) and U2{x,t), we prefer to write a self-contained equation, where the occasional 
amplifications follow from a purely stochastic dynamics that simulates the CML. More 
precisely, we introduce the model 

v{x,t) = {1 + Vl)w{x,t) , (9) 
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FIG. 1: Phase diagram of the RM model. In the two shaded regions ( / and //) the synchronized 
phase is unstable; in fact, above the solid line vp > 0. Along the dashed line, the TLE changes 
sign, so that region I corresponds to the linearly stable, but nonlinearly unstable regime, while II 
corresponds the linearly unstable region. Above A = Ac, the TLE vanishes together with vp- 



with 



w{x, t + 1) 



w 



.p. p = av{x, t) 



w(x, t + 1) 



av{x, t), w.p. 1 — p 



v{x,t)/A, w.p. p = aA 



if v{x,t) > A, 



(10) 



iiv{x,t) < A, 



av{x, t), w.p. 1 — p 

where a and A replace (1 — 2a)/ai and ^2/ (1 ~ 2cr), while w.p. is a shorthand notation for 
"with probability" . Only positive multipliers are assumed in order to guarantee a positive 
defined w{x,t) (simulations do confirm that the sign does not play a relevant role). Finally, 
periodic boundary conditions are assumed on a lattice of size L. 

The advantage of playing with this model is that it explicitly avoids the possibly subtle 
correlation that may be generated during the deterministic evolution of the CML and thereby 
spoiling the asymptotic behaviour of the observables we are interested in. Besides the 
probabilistic, rather than deterministic, choice of the amplification factor, the only other 
difference between the stochastic model ( p!ol) and the original set of two coupled CMLs is 
the distribution of the amplification factors that is dichotomic in the former case. We see 
no reason why this difference should affect the transition scenario. 

Moreover, in order to maximize propagation effects (that are responsible for the prop- 



agation of finite-size perturbations) we sliall restrict to the case e = 2/3 (the so called 
"democratic" coupling). Some rough numerical analyses do not, indeed, reveal qualitative 
changes when e is varied around 2/3. 

The most general way of testing the stability of the synchronized phase is by monitoring 
the evolution of a droplet of the unsynchronized phase. By denoting with N{t) the droplet 
size, i.e. the number of unsynchronized sites, at time t, the propagation velocity can be 
defined as 

. (11) 

i^oo 2t 

A negative TLE (the maximum Lyapunov exponent of model (p!0|)) implies that any in- 
finitesimal perturbation does decay. In spite of this linear stability, in Ref. it has been 



shown that vp can be positive, implying that the unsynchronized phase sustains itself and 
invades the synchronized one. By performing detailed simulations for different values of the 
parameters a and A, we have been able to construct the phase diagram plotted in Fig. 
The solid line, along which vp = 0, separates the synchronized from the unsynchronized 
phase (shaded region). The dashed line, along which the TLE is equal to 0, splits the un- 
synchronized phase into a linearly stable (J) and unstable (//) region. In the former one 
(ending approximately at A = Ac ~ 0.15), the nonlinear amplification mechanism prevails 
over the linear contraction induced by the negative TLE. Above Ac, the TLE changes sign 
exactly where vp vanishes too. 

Numerical analysis of stochastic synchronization in CML |^ suggests that when the 
TLE vanishes together with vp-, the critical properties of the synchronization transition are 
those of the MN class, while the transition is DP-like whenever vp only vanishes (the TLE 
remaining negative). 

Before entering a quantitative discussion about the nature of the transition in the present 
model, it is worth noticing a difference between the regimes I and II. Linear instability 
in II ensures that any finite perturbation of a synchronized state remains finite forever 
independently of the chain length. On the other hand, in J, a finite perturbation eventually 
dies in a finite chain. The reason why the synchronized regime can nevertheless be considered 
unstable is that the average life time of the perturbation diverges exponentially with the 
chain length. This is a typical property of systems in DP universality class and it can be 



traced back to the pecuhar nature of the "square root" noise amplitude in Eq. (p) [1 



A preliminary numerical analysis of the critical properties of the RM model for A = 




FIG. 2: Power law scaling behaviour in the RM model. In all graphs, the dashed lines correspond 
to the expected scaling behaviour (DP in a), b) and c); MN in d)). a) Absorption time as a 
function of system size for A = 0. Triangles corresponds to a = 0.6070, squares to a = 0.6063 and 
circles to a = 0.6051. b) Density of unsynchronized sites as a function of time for A = 0.1. The 
five solid lines correspond to (from top to bottom) a = 0.591, a = 0.59, a = 0.58955, a = 0.5893 
and a = 0.5890, respectively, c) Asymptotic density of unsynchronized sites as a function of the 
distance from criticality. Circles corresponds to A = 0.1 and squares to A = 0.2. d) Density of 
unsynchronized sites as a function of time for A = 0.2. The five solid lines correspond to (from 
top to bottom) a = 0.568, a = 0.5675, a = 0.5668, a = 0.5664 and a = 0.5662. All the graphs are 
plotted in a doubly logarithmic scales. 



and A = 0.01 has been already published in [|l^. Here we both perform more accurate 
simulations and extend the previous study to larger values (A = 0.1, 0.2) in order to find a 
signature of the change of critical behavior. In all cases, a is chosen to be the control pa- 
rameter, while the averaged (over different noise realizations) density pit) of unsynchronized 
sites will be the order parameter. The definition of p requires to fix a small threshold W to 
discriminate between synchronized {w{x,t) < W) and unsynchronized (w(x,t) > W) sites. 
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In principle, p{t) depends on W, both because the perturbation reaches different thresholds 
at different times and resurgencies can occur. Numerical analysis, however, indicates that, 
in practice, if W is chosen on the order or smaller than 10~^ no appreciable differences are 
observed. We shall come back to this problem in Sec. 4, to provide a more sound justification 
for the adopted procedure. 

In order to test the relationship between synchronization transition and the DP critical 
phenomenon, we have investigated the scaling behaviour in the vicinity of the transition. In 
DP it is known that, at criticality, the dependence of the density p{t) on t and L is described 
by the scaling relation [jl4| 



Pit) = L-'^g (J^j , (12) 

where z is the so-called dynamical exponent accounting for the dependence of the average 



time r needed for p to vanish with the system size L |19 



T ^ L' a = ac . (13) 

Since for small 9 = t/L^, the scaling function behaves as g{9) ~ 9^^, the exponent 6 turns 
out to describe the power-law decay of p{t) 

pit) r^r^ a = ac . (14) 

Finally, the exponent j3 characterizes the scaling behaviour of the saturated density of active 
sites po{t) as a function of the distance from the critical value, 

Po ~ (a — flc)^ a > ttc ■ (15) 

In analogy with usual nonequilibrium phase transitions, z, 6, and /? are expected to char- 
acterize all critical properties of the synchronization transition as well. In fact, simple 
dimensional arguments show that the exponents ruling the power law divergence exhibited 
by space- and time-correlation functions (while approaching the critical point) are linked to 
the previous ones by the standard relations 

z/|| = - , u^ = - . (16) 

Some of the scaling behaviors have been plotted in Fig. |^ to show the quality of the results, 
while a complete summary of the exponents are reported in Table |, together with the best 
known estimates for the DP and the MN |^ class. 



10 





A = 


A = 0.01 


A = 0.1 


A = 0.2 


DP 


MN 


z 


1.56 ±0.06 


1.58 ±0.02 


1.54 ±0.06 


1.5 * 


1.580745 ± 10"^ 


1.53 ±0.07 


5 


0.155 ±0.005 


0.15 ±0.01 


0.159 ±0.002 


1.2 ±0.1 


0.159464 ±6 - 10-6 


1.10 ±0.05 




0.24 ±0.02 


0.27 ±0.01 


0.27 ±0.01 


1.8 ±0.1 


0.276486 ± 6 • 10"^ 


1.70 ±0.05 




0.6063 . . . 


0.605. . . 


0.5895. . . 


0.5668 . . . 







TABLE I: Numerical results concerning the z the 6 and the /? exponent of the RM model are 
compared with the best available estimates. Values for ac indicate our best estimates of the critical 
point. Errors have been estimated as the maximum deviation from linearity in the log-log plot 
that it is used to extract the scaling law. * Indicates a value compatible with the theoretically 
predicted one, z = 1.5 (see text). 

The dynamical exponent has been estimated by averaging the behaviour of relatively 
small systems (from L = 2^ to L = 2^°) over a large number of noise realizations (of order 
10^). In order to minimize finite-size effects, the exponents 6 and (3 have been estimated 
from the time evolution of a single system of size L = 2^°, relying on the large size to reduce 
statistical fluctuations. In the MN context we have not been able to estimate z through 
the measure of the average synchronization time, but we verified, through finite size-scaling 



(|l^), that the value of the dynamical exponent is compatible with the theoretical prediction. 

Interestingly, similar results are obtained by adopting a different order parameter, i.e. the 
space averaged difference variable w(t) = {w{x,t))x. Also in this case, both {w(t)) (where (•) 
denotes ensemble average) and the absorption time ti, defined as the average time required 
for w{t) to become smaller than some threshold W, are found to follow the same critical 
scaling laws. The application of coarse-graining suggests that the space average is the 
"natural" order parameter in the context of both equilibrium and nonequilibrium critical 
transitions. 



III. FROM THE RM MODEL TO THE DP FIELD EQUATION 



In this section we investigate the connection between the RM model and the DP field 
equation (|). Let us first consider the simple case A = that corresponds to a discontinuous 
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but otherwise uniformly contracting local map. Eq. ([T0| ) can be recasted as 

w{x,t + 1) = 2av{x,t) - a^v^{x,t) + g{v)C{x,t) , (17) 
where ^'(x, t) is a zero-average 5-correlated noise with unit variance. In fact, 

e'(x,t) = ^[e.(a;,t)-(U] , (18) 



where 



and 



(a— l)f(x,t), w.p. 1— p, 



= (2a-l)t;-aV (20) 



If we now introduce the coarse grained variable p{x,t) = w{x,t) (where the bar denotes an 
average over a suitable space-time cell), we have that v{x,t) = p{x,t) + |V^p(x,t) so that 
Eq. yields, 

dtp{x, t) = aeV^p(x, t) + (2a — l)p(x, t) — a^p^(x, t) + a?ep{x, t)'V'^p{x, t) + 
<^[V^p(x,t)f + g{p{x,t) + IV'p{x,t))r^{x,t) . 



where, according to the central limit theorem |^2|], the coarse grained noise term ri{x,t) is 



Gaussian and 5-correlated in time and space. According to standard renormalizat ion-group 
arguments, [13], |14| the terms of order (V^p)^ and pV^p can be shown to be irrelevant. 



as well as the terms of order higher than or equal to p and a/V^p appearing in the noise 
amplitude g{p+ f V^p). 

From the definition (^1]) of g and after discarding the irrelevant terms, the above equation 
reduces to 



dtp{x, t) = aeV'^p{x, t) + (2a - l)p(x, t) - a?p^{x, t) + ^/ap{x,t)T]{x, t). (23) 

which is nothing but Eq. (^, thus proving that the synchronization transition in discontin- 
uous CML's can be traced back to a DP nonequilibrium phase transition. 

Let us now turn our attention to the more general case A > 0, which corresponds to a 
continuous local mapping. According to Eq. ([T0|), we have now to deal with two different 
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kinds of noise, depending whether v{x, t) > A or v{x,t) < A. By repeating the same formal 
derivation as in the previous case, we find that Eq. (|23| ) still holds when p > A, while for 
p{x, t) < A it must be replaced by the equation 

dtp{x,t) = — {2- aA)\/^p{x,t) + {2a-l- a'^A)p{x,t) + h{w)ri{x,t) , (24) 

where 

h{w) = p(x, t) - 3a2 + (2 + 5A)a3 - (3 + 2A) Aa^ + (1 - A2) Aa^ . (25) 

Accordingly, in Eq. (|^), the noise amplitude is proportional to the field itself, so that one 
should be led to the naive conclusion that the DP critical behaviour is destroyed as soon 
as A is finite, or, equivalently, that any CML system characterized by a continuous local 
mapping cannot exhibit a DP-like synchronization transition. However, the simulations 
described in the previous section suggest that DP-like transition can still been found for 
small but finite values of A. In the next sections we shall present theoretical arguments 
supporting such numerical findings. 



IV. FIRST PASSAGE TIMES 

In this section we clarify the problem of how and when it is possible to observe a DP-like 
scenario in models like the RM one, with no clearly identifiable absorbing state. As already 



noted in |T^, in any finite system (of length L) there always exists a finite probability for 
a generic configuration to be contracted forever, i.e. absorbed. A lower bound to such a 
probability is (in the discontinuous case), 

■ oo "I ^ 

Hil-WMa^ , (26) 

_n=l 

where 

wm = max w(x, 0) . (27) 

X 

However, since the null state, w{x,t) = 0, is reached in an infinite time, this configuration 
cannot be attained with perfect accuracy in numerical simulations and one is, in fact, obliged 
to fix a small but finite threshold. 
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The best way we have found to characterize the dependence of the perturbation evolution 
on its size is through an indicator closely related to the finite-size Lyapunov exponent (FSLE) 
introduced in Ref. |]23[. With reference to a perturbation initially set equal to 1 {w{x, 0) = 
1, X = 1...L), we introduce the first passage time Tg{W), defined as the (ensemble) 
average time required by the g-th norm of the state vector w(t), 



|w(<)IL 



i=l 



(28) 



to become smaller than some threshold W for the first time. At variance with Ref. |2^, 24 



we do not care if the evolution of the perturbation is non monotonous: as we shall see, in this 
context, the analysis does not only remain meaningful, but even more, it allows identifying 
the reason for the existence of a DP-like scenario even in the context of the continuous 
model. 

At variance with the standard Lyapunov exponent, the FSLE does depend on the choice 
of the norm (in particular, on the g- value in Eq. (|28|)). This circumstance is often considered 
as a difficulty, hindering a proper definition of FSLE: we prefer to see it as an indication of 
a richer class of phenomena associated with the evolution of finite-amplitude perturbations. 
Anyway, it has been noticed in Sec. || that the "natural" order parameter of DP transition is 
the spatial average of the state vector. Accordingly, we have decided, in the present context, 
to fix g = 1 (that corresponds to performing an arithmetic average) and to drop, for the 
sake of simplicity, the dependence on q. 

The FSLE A(W^) can be introduced by first fixing a sequence of decreasing thresholds 
Wn,n = 0,1, 2,..., 

= r , r < 1 , (29) 



n-l 



and by then defining 



Inr 



T{Wn+l) - T{Wn) 

where the dependence on the "discretization" Inr is left implicit. In the limit r 
definition becomes 

A(W) ''^<'*'' 



-1 



(30) 

1, the 

(31) 



_d{\nW)_ 

In the further limit W 0, A{W) reduces to the usual Lyapunov exponent A, independently 
of the adopted g- value. When A = 0, A = Ina. 
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A. Uncoupled limit 



In the uncoupled limit, e = 0, each site converges independently to the synchronized state 
(as long as a < 1). In spite of the low- dimensionality of the problem, even in this case, an 
analytic expression for the FSLE can be obtained only at the expense of introducing further 
approximations. We shall see that the resulting expression can be nevertheless profitably 
used even in the coupled regime. 

By setting r = a and restricting ourselves to the case A = 0, it easy to show (see appendix 
0) that 

where t{Wq) = 0, Wn = a"'. By inserting Eq. ( P^ in Eq. (^) one obtains 

Eq. (|32D implies that, for ^ oo, 

r(W„) ~n + no = i^+no . (34) 
ma 

By inserting Eq. (|3^ ) into Eq. (|33D and recalling that A = Ina (for A = 0), we obtain 

^(^") = 1 + anoWr. + (a/lVn In l^n^ ' ^^^^ 
As we are interested in describing the region where Wn -C 1, and owing to relative smallness 
of a (a <^ uq), this equation can be further simplified to 

^<'^' = rrw^w^TF ■ 

where we have also dropped the unnecessary dependence on the index n. In the W ^ 
limit, the leading correction to the standard Lyapunov exponent is provided by the term 
proportional to bi. From the structure of Eq. (|5BD, it is natural to interpret the inverse 
of bi as the critical threshold Wc, below which the dynamics of the uncoupled system is 
dominated by the maximum Lyapunov exponent A. 

From Eq. ( p6D and by integrating Eq. (^T]), it is possible to derive also an analytic ex- 
pression for the first passage time t{W), 



t{W) 



j AiW)~^d{\nW)^j[lnW + {bo + bi)W-biW\nW] + b2 , (37) 



15 



where 62 is the integration constant. In principle, one could imagine of determining 62 by 
imposing r(l) = 62 — bi/X equal to 0, since the evolution starts precisely from W = 1. 
However, we cannot expect our perturbative formula to describe correctly the initial part of 
the contraction process. Therefore, 62 must be determined independently. 



B. General case and scaling arguments 

In the coupled case, we have not been able to derive an analytical expression for the 



FSLE. Nevertheless, a comparison with the numerical results has revealed that Eqs. (p6| , |37D 
describe in a reasonable way the dependence of A and t owW even in the continuous model. 
However, while A still denotes the standard Lyapunov exponent of the process and can 
thus be computed independently, now 60, 61, and 62 have to be determined by fitting the 
numerical data. We have preferred to keep also the term proportional to &o (relevant only 
for relatively large VT- values), since its presence guarantees a much better reproduction of 
the numerical data. In Fig. we see that Eq. ( |37D provides a good parametrization of the 
numerically determined r-values over a wide range of thresholds, both in the discontinuous 
and the continuous model (see the solid curves). In panel b) we notice that, although the 
theoretical expression (|36|) does not provide an equally good description of the FSLE, it is 
nevertheless able to pinpoint the crossover towards the small-amplitude behaviour of the 
perturbation. We will see that the possibility of identifing the largest scale (defined by I/&1) 
over which the linearized dynamics (described by the standard Lyapunov exponent) sets in 
represents a crucial point of our analysis. 

It is now natural to ask to what extent Eq. ( pTl) is able to account for the scaling behaviour 
in the vicinity of the transition. By replacing p with W in Eq. (|T2|) and t with the first passage 
time r, one expects that, at criticality, 

W = L-^'g{t/V) . (38) 

Inversion of this equation leads to 

t{W) = L'g-\WL^') . (39) 



Before mutually comparing the two expressions (|37| , |39|) , it should be first stressed that 



they have been introduced to address different questions. On the one hand, Eq. ( ^7|) is 
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FIG. 3: Numerical data for the first passage time (left panel) and the FSLE (right panel) at the 
critical point are fitted with Eq. ( |37| ) and Eq. (^) (full lines). Circles refers to A = 0.01, L = 256, 
flc = 0.6055 and squares to A = 0, L = 128, Oc = 0.6063. 

an approximate expression introduced to account for the crossover towards the VT-range 
where the dynamics is controlled by linear mechanisms and no scaling behavioiur should be 
expected. On the other hand, Eq. (^) is a rigorous but implicit statement about the scaling 
region only. 



Compatibility between Eqs. (pTf) and ( pOf ) requires a proper dependence of 6o, and 62 
on the systems size L, namely 



ho = -A 
62 = W , 



60 - ^1(1 + 6zhiL) 



(40) 
(41) 
(42) 



where Bq, bi, and 62 are suitable positive constants. By inserting Eq. (^) into Eq. (0), one 
finds that 



Inl^ 

r{W) = ^-L^ 
A 



boWL'^ - biWU^ \n{WL^') + 62 



(43) 



from which we see that the first term in the r.h.s. is the only one which does not follow the 



required scaling law (|39D . In fact, (lniy)/A, accounts for the linearly stable behaviour in 
a regime where a finite-state model (such as, e.g., the famous Domany-Kinzel model [jlBp 
would be otherwise characterized by a perfect absorption (when a configuration of all O's is 
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FIG. 4: Finite-size scaling behaviour at the critical point of both h\ and 62 for A = and = 
0.6063 (left panel) and A = 0.01 and Oc = 0.6055 (right panel). The dashed lines indicate the 
best power law fit. Left: 61 (circles) scales as L^''^^, while 62 (triangles) scales as L^'^^. Right: bi 
(squares) scales as L^-^'^, while 62 (diamonds) scales as L^-^^. 

attained). 

In order to test the correctness of the whole picture, we have studied the dependence of bi 
and 62 on L. In Fig. their behaviour is plotted at criticality for the discontinuous and the 
continuous model: both quantities show a good agreement with the power law divergence 
predicted by Eq. ( ^2]) {z ^ 1.58 and z{l + 6) ~ 1.82). As for the last parameter bo, given 



its involved dependence on L and the approximate character of Eq. (43), we can only claim 
that its dependence is quahtatively consistent with the theoretical prediction. 

One of the most important results of our study is the objective identification of a threshold 
Wc = 1/1^1 1, below which linear stability analysis holds and its scaling dependence on L 
[Wc ~ L~^). In a model like the cellular automaton considered by Domany and Kinzel [IB 



absorption in a finite system occurs when all sites collapse onto the absorbing state: this 
means that the minimal meaningful density that can be considered is pm = 1/L. In the 
present context, Wc plays the role of pm'- below Wc, the critical behaviour is dominated 
by the linearly stable dynamics. The difference between the two systems lies in the scaling 
dependence of the maximal resolution on L. Since Wc decreases faster than 1/L this means 
that, e.g., the scaling range for Wc is wider in the present model than in finite-state systems. 
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Finally, we comment about the reason why the range of validity of the linear stability 
analysis can eventually vanish even in models like the continuous RM, where every pertur- 
bation locally smaller than A should behave linearly. The reason is that t{W) is defined 
as the average first-passage time: even if the perturbation is homogeneously small, if L is 
sufficiently large, some occasional amplification may occur and drive, on the average, the 
system out of the linear region. It is only below Wc that such sporadic resurgencies are 
sufficiently rare not to modify significantly the stable linear behaviour. 



V. CONCLUSIONS AND OPEN PROBLEMS 



In this paper we have expounded a partially rigorous argument to show why the synchro- 
nization transition in spatially extended systems may belong to the DP universality class. 
Although our theoretical considerations restrict to the discontinuous RM model, scaling 
analysis of the first-passage time t{W) suggests that the transition belongs to the DP class 
also in a finite parameter region of the continuous model. Since direct numerical simulations 
in the more physical class of CMLs have been basically restricted to discontinuous maps, we 
find it wise to test the validity of our conclusions also in the context of continuous, though 
highly-nonlinear maps. Accordingly, we have considered two lattices of maps coupled as in 
Eq. the local map is chosen similar to those defined by Eq. (0), namely 

x/ai < X < ai 

1 — (x — ai){l — a3)/a2 ai < x < ai + a2 ■ (44) 
tts + a4{x — a\ — a-i) ai + 0:2 < x < 1 

with «! = 1/2.7, 0:3 = 0.07, = 0.1. The reason for this choice is that in Ref. it 
has been shown that in such a model (for the same parameter values and 0:2 < 0.013 p6[ ) 
nonlinear effects prevail over linear ones. In fact, it was observed that the propagation 
velocity vp of finite- amphtude perturbations (see Eq. ([TTl) ) is larger than the propagation 
velocity of infinitesimal perturbation (for a definition of vl, see ^7, ^). For instance. 



fix) 



for 02 = 4 ■ 10~^ and e = 2/3, vl = 0.4184, while vp = 0.5805. In this regime, upon 
varying the coupling strength a, synchronization arises through a continuous phase transition 
accompanied by a negative transverse Lyapunov exponent and a vanishing vp d.t the critical 
point (Tc = 0.17756 .... As it can be appreciated in Fig. ||, where we have plotted the space 
averaged difference variable w{t) versus time for different values of the control parameter, 
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FIG. 5: Log- log plot of the space averaged difference variable w{t) as a function of time for two 
coupled CML's (see text) and for different coupling values: from lower to upper, the full lines 
correspond to cr = 0.178, a = 0.17756, a = 0.1775, a = 0.177, while the long dashed line marks the 
power low expected for the DP critical behaviour. Numerical data have been obtained averaging 
over 100 realizations of a CML of size L = '2P . 

the critical decay rate is 5 = 0.158 ± 0.01, fully compatible with the expectation for a DP 
transition. We are thus reinforced in the conjecture that the DP scenario is robust and not 
just restricted to the highly nongeneric case of discontinuous maps. 

The crucial difficulty to determine the universality class for the synchronization transition 
is that the order parameter (the difference field) can be arbitrarily small. This casts doubts 
about the very definition of the zero-difference field as a truly absorbing state. In fact, in 



a previous paper ||T5[, it was speculated that the DP scaling behaviour might be restricted 
to a finite range. The analysis carried on in this paper clarifies that the synchronization 
transition genuinely belongs to DP universality class: this has been understood from an 
objective identification of the threshold below which the dynamics is really controlled by 
linear mechanisms and thus corresponds to an effective contraction. The parametrization 
of r(I^) introduced to describe the single-map case has greatly helped to unveil the overall 
scenario since it has clarified that the basic effect of the diffusive coupling is to renormalize 
the parameters defining TiyV) (see Eq. (pT])). Here, the parameter values (in particular VFc) 
have been inferred by fitting the numerical data; in the future, it will be desirable to find 
an analytic, though approximate, way of performing the renormalization. 

Once we have concluded that synchronization arises through a DP-like transition in a 
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finite parameter region, it is natural to ask how this scenario crosses over to the standard 
transition characterized by a vanishing of the Lyapunov exponent and by the KPZ critical 
exponents. With reference to Fig. |I], this question amounts to investigating the region around 
the multicritical point A^. A purely numerical analysis of this region is not feasible in this 
model, as it would require considering too large systems to be effectively handled. We are 
currently studying this problem in a different context, where preliminary studies indicate 
the possibility to draw quantitative conclusions. 

Finally, since it is known that finite-size Lyapunov exponents do depend on the norm, it 
might be worth considering q values different from 1, in order to check to what extent the 
universality of the transition is preserved when different averaging procedures are adopted to 
assess the amplitude of the global perturbation. In particular, since q = oo (corresponding 
to the maximum norm) takes care only of the extreme fluctuations of a perturbation field, 
it is not totally obvious that the behaviour of the corresponding first passage time follows 
exactly the above described scenario. 
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APPENDIX A: FIRST PASSAGE TIMES IN THE DISCONTINUOUS UNCOU- 
PLED RM MODEL 

In this appendix we report the analytical calculation of the first passage time when e = 
and A = 0, to prove Eq. (^2]). Being Wn = a" and w{0) = 1, we have also t{Wo) = 0. In 
order to compute the first passage time through a threshold Wn, we need to know the average 
time needed to pass from Wn-i to Wn- With a probability 1 — aWn-i, this can occur in 
one time step, if the amplification mechanism is not activated and the synchronization error 
is contracted by a factor a. On the other hand, with probability aWn-i, the amplification 
resets the state variable to the value 1. In this case, one has to wait for the synchronization 
error to shrink back to the n-th threshold, which, by definition, occurs in an average time 
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Tn-i- At this point, the error can either shrink to Wn or be reset again to 1, to start again 
the process. Altogether, 

T{Wn) = T{Wn-l) + 1 • (1 - aWn-i) + (2 + t(W^„_i)(1 - aWn-i)aWn-i + 

(3 + 2T(Wn-l)(l - aWn-i)(aWn-i)'' + ... 

= r(P^„_i) + (1 - aWn-i) E,=o [(1 + ^ + tr{Wn-i))iaW^-iy] 
= riWr^-i) + (1 - a-) [E.=o(«")^ + iriWn-i + 1) E.=o ^ ■ (Al) 

Summing up the series, one obtains 
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